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Q ■ Abstract 

We study the rational potentials V(x), with sextic growth at infin- 
ity, such that the corresponding one-dimensional Schrodinger equation 
has no monodromy in the complex domain for all values of the spec- 

^ ■ tral parameter. We investigate in detail the subclass of such potentials 

which can be constructed by the Darboux transformations from the 
well-known class of quasi-exactly solvable potentials 

6 2 Kl + 1) 

V = x° - vx l + 



X 2 

We show that, in contrast with the case of quadratic growth, there are 
monodromy-free potentials which have quasi-rational eigenfunctions, 
but which can not be given by this construction. We discuss the rela- 
tions between the corresponding algebraic varieties, and present some 
elementary solutions of the Calogero-Moser problem in the external 
field with sextic potential. 
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1 Introduction 

Consider the Schrodinger equation with a meromorphic potential V(x): 

(-D 2 + V(x))ij = Xij, (1) 

where here and below D — j-. We say that such a potential V(x) has 
trivial monodromy, if all the solutions ip of the corresponding equation (TjQ) 
are meromorphic in x for all A. The general problem is to describe all such 
potentials. 

In such generality this is probably a hopeless question. However if we 
restrict ourselves to a suitable class of potentials it is known to be solvable. 
The first results in this direction were found by Duistermaat and Grunbaum 
[lj, who solved it in the class of rational potentials decaying at infinity: 

_ y^ kjjkj + 1) 

^ (x-Xi) 2 ' 

They showed that the corresponding parameters ki must be integers, and 
all such potentials are the results of Darboux transformations applied to 
the zero potential. Therefore they are given by the Burchnall-Chaundy (or 
Adler-Moser) explicit formulas [21 [3j. The corresponding configurations of 
the poles Xi are very special: in the case when all the parameters k t = 1, 
they are nothing else but the (complex) equilibria of the Calogero-Moser 
system, with the Hamiltonian 

N N 



Jb j t\j n 



which are described by the following algebraic system: 

N 

[Xi Xj 



- 1 

A remarkable fact (discovered first by Airault, McKean and Moser [1]) is 
that this system has no solutions unless 

m(jn + 1) 



is a triangular number, in which case the solutions depend on m arbitrary 
complex parameters. 



Oblomkov [5] generalised the Duistermaat-Griinbaum result to the case 
of rational potentials with quadratic growth at infinity: 

y _ x 2 + y^ k i{kj + !) 

^-f (X - Xi) 2 ' 
1=1 v ' 

by proving that all such potentials can be found by applying successive Dar- 
boux transformations to the harmonic oscillator, and that they have the 
following explicit description in terms of the Wronskians of Hermite polyno- 
mials Hk(x): 

V = x 2 - 2D 2 \ogW(H kl ,H k2 , . . . , H kn ). 

In particular the corresponding Schrodinger equations must have infinitely 
many solutions of the quasi-rational form 

t/j = R(x)exp{±x 2 /2), (3) 

where R(x) is some rational function. 

An interpretation of this result in terms of the generalised Stieltjes rela- 
tions and the rational solutions of the dressing chain [7] was proposed in the 
paper [6] , where also a new class of meromorphic monodromy-free potentials, 
expressible in terms of the Painleve transcendents, was found. 

In this paper we consider the set M of monodromy-free even rational 
potentials with sextic growth at infinity: 



" = «(*) + E^^r W 



1=1 



where P&(x) is an even polynomial of degree 6, which for simplicity will be 
assumed to be of the form 

P 6 (x) =x 6 - ux 2 , (5) 

where v is a parameter. 

Our motivation came from the theory of Calogero-Moser systems. It is 
known, after Inozemtsev [8], that the corresponding classical system in an 
external field with the Hamiltonian: 

N N AT 

, n ^ (Xi- Xj) 2 

1 = 1 2=0 ]^l V J 

is integrable (note that this is believed to be not the case if the degree of 
P is higher than 6), so one can expect some effective way to describe the 
equilibria of this system and the corresponding set M. 



Having in mind the case of quadratic growth, let us consider also the set 
S of potentials (jlj), such that the corresponding Schrodinger operator has a 
quasi-rational eigenf unction, of the form 

ip = R(x)exp{±x i /£), (6) 

where R(x) is some rational function. 

Finally we consider the set D of the operators, which are Darboux-related 
to the operators of the form 

o k n 9 mini + 1) 

H P , m = -D 2 + x 6 - (2p + 1 )x 2 + K } (7) 

x z 

with integer p and m. The operators £ Pjm are known, after Turbiner, P, [TO] , 
to be quasi- exactly solvable, and they were studied later from various points 
of view in [TTJ [T2J [13], [TU [151 HE]- They have M quasi-rational eigenf unctions 

# m) = R { t m \x) exp(±s 4 /4), % = 1, . . . , M, 

where M is a positive integer determined by p and m (see section 2 below). 
The set D consists of the potentials of the form: 

V(x) =x Q - (2p + l)x 2 + m(m 9 + X) - 2D 2 log Wj(x), (8) 

x l 

where I C {1,2, ...,M} and Wj is the Wronskian of the corresponding 
quasi-rational eigenfunctions fy, i € I. 

We will study all these three sets and the corresponding affine algebraic 
varieties. In particular, we show that, under some assumptions, 



"■) 



D C S C M, 

and that, in contrast to the quadratic case, all the inclusions are proper: 

D + S + M, 

even if we restrict the definition of the set M to odd integer values of the 
parameter v. This means that only a proper subset of monodromy-free po- 
tentials with sextic growth at infinity have an explicit description in terms 
of Darboux transformations. 

We study also the corresponding time-dependent Schrodinger equation, 
and show that the Stieltjes relations for them can be considered as Backlund 
transformations between Calogero-Moser systems in different sextic poten- 
tials, with different numbers of particles. We construct a family of explicit 
solutions in elementary functions for such a system. We note that these solu- 
tions are very special - Inozemtsev's general solution of the system is highly 
transcendental, being expressed in terms of Riemann ^-functions. 



2 The locus conditions 

In this section we recall the trivial monodromy conditions pQ, which we will 
call locus equations following the terminology of [I]. 

Suppose that a Schrodinger equation with a rational potential V(x) 

(-D 2 + V(x))ij = X^j, 

has all its solutions meromorphic for all A. At a singular point, which we 
may take to be x = 0, we have the two series: 

oo 

V(x) =X- D J2 V n% n , 
n=0 
oo 

^(x,\)=X- m Y,*PnX n . 



n=0 



where we suppose D,m > 0, and ip ^ 0. Then the Schrodinger equation 
becomes: 

oo oo 

n=0 1=0 

oo 

= Xx~ m J2 ^nX n . 
n=0 

Expanding, we find: 

_. oo oo n 

-- Y, (n + 2)(n + l)^ n+2 x n + x^ m+D ^ xn Y,^ V ^ 

n=-(m+2) n=0 1=0 

oo 

= Xx- m J2^nX n . 
n=0 

Hence, to balance the leading terms, it follows that V(x) has only double 
poles, D = 2, and that V-2 — \rn{m + 1). The integer m will be called the 
multiplicity of the pole. 

There are further constraints. The coefficient of x n includes the term 

-{m{m + 1) — n{n — l))^ n 

and other terms which involve ipiV n -i, for I < n. Hence the m+ 1st equation 
does not involve tp m+1 . Thus the Schrodinger equation can be solved to all 



orders, provided that the (m + l)st equation is consistent with the previous 
equations. The resulting system may be written in matrix form: 



( ° 



























2m -1 


m 


v 1 


Ui 


v 2 -\ 



m V\ \ 

2m -I V t V 2 -X 

ml- 1(1 -l)/2 Vi ... Vi 
Vi V 2 - A ... 



Vo-\ 



V2m-1 

V 2m 
V 2m V2m+1 / 



( i>2m \ 
■02m- 1 



^1 
V ^0 / 



/o\ 





w 



The (2m + 1) x (2m + 1) matrix on the left must be singular for all values of 
A; it follows from this that the leading m + 1 odd coefficients of the Laurent 
series of V all vanish: 



V 1 = V 3 = ... = V 2m+1 = 0. 
In the special case of potentials of the form 






(9) 



>iV 



where P(x) is some polynomial, we have an algebraic system of K = ^ 3 =i kj 
equations in iV variables xf. 



N 



kj{kj + 1) 



P^Hx,) - (2s)\J2 ( AJ \ 2+1 = 0, i = l,...,N,s = l,...k i . (10) 
If all the multiplicities hi — 1 we have N equations for N variables: 

N 



p'ixi) - E 



3+i 



{Xi Xj ) 



A3 



0, i = l,...,N. 



'ir 



In our case, with P(x) = x 6 — ux 2 , we will consider another special 
case, when x = may have an arbitrary multiplicity I, all other poles have 
multiplicity 1, and the potential is even: 



N 



V(,r)^,r"-z/x 2 + ^±il + E 



x/- 



^ (x - Xif 



(12) 



where xi, . . . ,xn, N = 2n is a set of distinct non-zero complex numbers 
symmetric with respect to the origin. In that case the locus equations at 



x = are automatically satisfied, and we have the following algebraic system 
on xf 

v^ 2 1(1 + 1) „ 

^ J/ * 

This is an affine algebraic variety NL^ N '(v, I), which is non-empty for all values 
of N and the parameter v. 

What we are going to show here is that, for special values of the parameter 
v, there is a subset of such equilibria which can be described quite explic- 
itly. They correspond to a special subclass of the monodromy-free potentials 
(HI), and form some zero- dimensional algebraic varieties, whose properties we 
would like to study 



3 The Darboux transformation 

We recall here the classical procedure going back to Darboux |17| . 

Let us consider the Schrodinger operator L^ with potential V^(x), its 
eigenf unctions ip(°\ and in particular its eigenfunction ip[ with eigenvalue 
Ai: 

(L (0) - Ai)V (0) = {-D 2 + V (0 \x) - Ai)V> (0) = (A - Ai)^ (0) - (14) 

The Schrodinger operator can be factorised as 

(-D 2 + V {0 \x) - Ai) = -(D + v^)(D - v^). 

The operator on the right must annihilate the eigenfunction ift[ , so we may 
take i/ -* = D\n.(il)\ ). The potential V^(x) is given by 

V (o) {x) _ Xl = ^l + v (o) 2 . 
ax 

Now we act on the left of (IT3I) with the operator (D — v^). This takes any 
eigenfunction ip^ of L^ into a new function 

^W = (D - t,(°>)^°) 

which is either zero, if ip^ is a multiple of ij)\ , or else an eigenfunction of a 
new operator L^\ Explicitly 

-{D - v^)(D + v^)(D - i/°>)^°> = (A - X )(D - v^)^°\ 



so that 

(L (1) - Ax)^ (1) = -(D - v<®)(D + v {0) )^ = (A - Ax)^ (1) . (15) 
The potential V^(x) in the transformed operator 

L<V-\ 1 = -{D-v<- ')){D + v<®) 
is then given by 

ax 
Thus 

VM(x) - V<®{x) = -^— = ~hi4 0) . (16) 

ax ax 2 

If we now choose one of the remaining eigenvalues A2, say, then the trans- 
formed Schrodinger equation, shifted by this eigenvalue, is 

(L (1) - A0^ (1) = (-D 2 + V {l \x) - Ai)^ (1) = (A - Ai)'0 (1) . (17) 

We may now repeat this procedure indefinitely. The corresponding operators 
L^ are related by 

where the n-th order 'dressing operator' is given by: 

£>W = {D- v {n ~ l) ){D - v {n ~ 2) ) ...{D- v w ). 
Since D^ must annihilate the chosen set of eigenfunctions 

of L(°\ with respective eigenvalues Ai, A 2 . . . , A n , we have 



^(n) = D (n)^(0) 



wM \4 \...^ ] ) 



where W n denotes the Wronskian of its n arguments. Further, it may be 
shown by induction (see Crum |18j ) that 

V V>(x) = V^(x) - 2^1n W^i 0) , V4 0) , • • • ,^ 0) )- (18) 

What is important for us, is that if the initial potential V^ is meromor- 
phic and has trivial monodromy, then it follows from fll8p that the same is 
true for V^ as well. We will use this fact in the next section to describe a 
class of monodromy-free potentials with sextic growth at infinity. 
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4 The even sextic potential 

The anharmonic oscillator in an even sextic potential, with Hamiltonian 

H= hp 2 + x 6 -vx 2 ) 
has the special property that the corresponding system 

d2x a 5^ 

— — = — OX + vx 

at 2 
can be solved explicitly using the elliptic Weierstrass p-functions: 

1 



x(t) 



p(V2E(t-to);g 2 ig 3 )-7k 



where the invariants g 2 and g% are given by: 

v 2 

92 ~3E^ 
2 z/ 3 

93 ~E~ 27£ 3 ' 

and the energy E and starting time to are the two arbitrary constants. 

In contrast, the Schrodinger equation for the corresponding quantum sys- 
tem 

(-D 2 + x 6 - vx 2 ) if> = Xip (19) 

can not in general be solved exactly; instead it is at best quasi-exactly solvable 
in the sense of Turbiner [TU] . 

To be more precise, let us consider the following slightly more general 
family of Schrodinger equations pH] 



(.£2 +x 6_ ux 2 + l l±R)^ = ^ (20) 

X 2 

studied in detail in particular in (TTj [TJ1 Ell HH EE1 HE] • Let us look for a 
special class of quasi-polynomial solutions of the form: 

ip = x"P(x) exp(±s 4 /4), (21) 

where P(x) is a polynomial. Let us consider first the case when 

ip = x l "P(x) exp(-x 4 /4). (22) 

9 



It is useful to rewrite this equation in terms of 0, where ip = <p(x) exp(— x 4 /4): 
-(D 2 - 2x 3 D - (V - 3)x 2 + ^4^)0 = A< ^ 

JL 

The coefficients a n of a series expansion for <p\ 

oo 

(j) = X^Y^ a nX U , 
n=0 

satisfy the following recurrence relation 

[(^ + n + 2)(fi + n + l)-l(l + l)]a n+2 + \a n + (u-2n-2fi + l)a n _2 = 0. (23) 

The very first equation gives 

fi(jjt - 1) - 1(1 + 1) = 0, 

which implies that either /i = / + 1 or jjl = — I. In the first case, the recurrence 
becomes: 

(n + 2) (n + 21 + 2)a n+2 + Xa n + {v-2n-2l- l)a n _ 2 = 0, (24) 

in the second: 

(n + 2){n-2l + l)a n+2 + \a n + {v-2n + 2l + l)a n _ 2 = 0, (25) 

In particular, we see that the even and odd-index terms are decoupled, 
which is a consequence of the symmetry of the original problem. The quasi- 
polynomial solutions exist when v — 2n — 2/ — 1 = for some even positive 
integer n: in that case the a n can all be set to zero for n > 2M, where 

= ,-(21 + 1) 



which in that case is a positive integer. The eigenvalues A must satisfy the 
characteristic equation of the corresponding tridiagonal matrix, which is an 
algebraic equation (with integer coefficients if/ is an integer) of degree 2M—2. 
These polynomials are studied in more detail by Bender and Dunne [H] and 
sometimes called Bender-Dunne polynomials. 

When /i = — / we have M quasi-polynomial solutions (J22l) . where 

M= " + (2 ' + 1) 
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provided M is a positive integer. 

Similarly, the quasi-polynomial solutions of the form 

^ = x l+l P(x) exp(a; 4 /4), (26) 

exist when 

M= -»-( 2l + l) 
4 

is a positive integer; in that case there are M such solutions. For the quasi- 

polynomal solutions of the form 

= aT'P(x)exp(x 4 /4), (27) 

the same is true for 

- v +(2 + 1) 



4 
Summarising the above, we have: 

Proposition 4.1 The Schrodinger equation $20\) with integer I may have 
quasi-polynomial solutions of the form (21\) only if v is an odd integer. In 
that case the number of such solutions is M, which is the non-negative integer 
of the form 

±v ± (21 + 1) 

M = \ ■ 28 

4 

The combination of the signs in this formula for M determines the type of 
the corresponding solutions l[21\) . where P(x) must be an even polynomial of 
degree 1M — 2. The corresponding eigenvalues Ai, . . . , Xm are distinct real 
algebraic numbers. 

In particular, we see that such solutions exist for all odd v except v = 
±(2/ + l). 

Example. Let I — 0, v — 7: 

L = -D 2 + x 6 - 7x 2 . (29) 

Then M = ^±1 = 2. Both signs are plus, so the corresponding quasi- 
polynomial solutions are of the form 

■0 = x~ l P{x) exp(— x 4 /4) = P(x) exp(— x 4 /4), P(x) = a + a 2 x 2 . 

The characteristic equation is 



A 2 
4 A 



0, (30) 
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giving A = ±2a/2. For A = — 2v^2, the eigenfunction is found to be 

i> x = (V2x 2 + 1) exp(-x 4 /4) (31) 

which, having no zeroes, must be the ground state, while for A = 2v / 2, the 
eigenfunction is 

ijj 2 = (V2x 2 - 1) exp(-x 4 /4). (32) 

This has two zeroes, so there must be another (odd) excited state, with a 
single zero, with eigenvalue between in the interval (— 2a/2, 2y / 2) - however 
this eigenfunction is not quasi-polynomial. 

5 Dressing the sextic potential 

The potentials 

6 2 ^+!) 

V = x b - vx 1 + 



.i' 2 



with integer I are clearly monodromy-free, so we can apply the Darboux 
transformations to construct more operators with trivial monodromy. If we 
use only quasi-polynomial eigenfunctions, all the new potentials will be ra- 
tional with sextic growth at infinity Indeed by the Crum formula fll8p the 
new potential is 

V I = V-2D 2 lnW I (x), 

where Wi = W(ipi 1 ,ipi 2 i ■ ■ ■ > V'vJ is the Wronskian of the quasi-polynomial 
eigenfunctions (I2T1) from a subset I — {ii,i 2 , ■ ■ ■ , i m } C {1, . . . , M}. 
In the case when these eigenfunctions have the form 

ipi(x) = x l+1 Pi(x) exp(x 4 /4), (33) 

we have the following 

Proposition 5.1 The Wronskian Wi has a form 

Wi = a , m(i+1)+ i m(m _i )p ^^ exp ( mx 4/ 4 ) ( 34 ) 

with an even polynomial Pi of degree 

2d = 2m(M - m) (35) 

non-vanishing at zero. 
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Indeed we have 

Wtyi^, • • • , VO = x m{l+l) exp(mx 4 /4)^(P il , P i2 , . . . , P lm ), 

where Pi are the corresponding even polynomials whose coefficients a^ are 
determined recurrently by (1231) . We claim that 



W(P n ,P l2 , ...,P lm )= x^ m -Vp(x) (36) 

for some even polynomial P with -P(O) 7^ 0. To prove this let us assume with- 
out loss of generality that ik = k. We can apply elementary row operations to 
the coefficients of the polynomials Pi, ... , P m so that the lowest order terms 
are of increasing degree: 

P k = x 2k Q k (x) 

for some even polynomials Qk(x). We claim that Q k (0) 7^ 0. For this we need 
to prove the independence of the corresponding vectors 

a = \ a > a 2 5 • • • 5 a 2/J> % = 1; • • • ) ^ 



determined by the recurrence relations (1241) . From these relations it follows 
that <i2j is a polynomial in A of degree j and thus a^j are the values of these 
polynomials at the eigenvalues Aj. Since we know that all Aj are distinct the 
claim now follows from the Vandermonde formula. 

This implies that the leading term of the Wronskian 

W(P il ,...,P im ) = W(P il ,...,P im ) 

at zero is 

0+l+...+(m-l) _ x \m(m-l) 

Applying a similar procedure at infinity we can instead make the degrees 
strictly decreasing: 

degPi = 2M - 2«, 

so the degree of the Wronskian is 

2M-2 + (2M-5) + (2M-8) + . . . + (2M-3m+l) = m(2M-2) — m(m-l). 

This means that the degree of the polynomial P(x) in ( |36j) is 

3 1 

m(2M — 2) m(m — 1) m(m — 1) = 2m(M — m). 

This completes the proof. 

13 



Thus the result of the corresponding Darboux transformations has the 
form 

R / x 9 1(1 + 1) + 2m(l + 1) + m(m - 1) „ 9l „ , . 

x 2 

which is equal to 

« s 9 U + m)(l + m + 1) 9 , . . . , 

Vj = x 6 - (u - 6m)x 2 + K JK 7 —^-2D 2 \nP I (x) 37 

with an even polynomial Pj of degree 2d = 2m(M — m), where 

M = "- (2l + 1) 

4 

We will denote this set of potentials -D+- having in mind the combination of 
signs in the formula for M. 

Similarly we have three more sets: 

D ++ consisting of the results of m Darboux transformations using the 
eigenf unctions of the form ipi(x) = x~ l P(x) exp — x A /A : 

Vl = x 6 - (v - Qm)x 2 + ( ^~ m)(/ 2 ~ m + 1) - 2D 2 \nP I (x), 



degPi{x) = 2m(M - m), M 



x z 

u + (2l + l) 



4 

D |. consisting of the results of m Darboux transformations using the 

eigenf unctions of the form ipi(x) = x~ l P(x) expx 4 /4 : 

Vj = x 6 - (u + Qm)x 2 + (/ ~ m)( ^~ m + 1) - 2D 2 InP^x), 

x 1 

degPx(x) = 2m(M -m), M = ~^ . + ; 

D consisting of the results of m Darboux transformations using the 

eigenf unctions of the form ipi(x) = x l+l P(x) expx 4 /4 : 

V I = x 6 -(u + 6m)x 2 + ( l + m )( l + m + 1 ) _ 2D 2 lnP/(;E ) 5 

x 2 

degPx(x) = 2m(M -m),M= ~ U ~ ^ — iLJ.. 

However there are two natural duality relations between these sets. 
They correspond to the spectral equivalences between the operators ( 1201) 
studied by Dorey, Dunning and Tateo in [16J. 

14 



First of all the form of the potentials in (120]) is clearly invariant under the 
change 

l->-l-l,u-+u, (38) 

which means that the sets 

D ++ = D + _, D_ + = D__ 

provided we allow / to be any integer. 

The second duality is more interesting. Let us note that in the D + _ case, 
when m = M, the resulting potential is 

V M = t/W-(„-6A/)^ + ( ' + M)( ' + M+1) (39) 

X 2 

which corresponds to the change 

^(/* = (/-6M, l^l* = l + M. (40) 

We can formally add 

M -f M* = -M (41) 

to make it an involution, which is nothing else but the "third spectral equiv- 
alence" from |16j . 

Note that we can allow M to be a negative integer, provided we replace 
M by \M\ in the formula for the degree of the polynomials Pi 

2d = m(\M\ -to). 

This duality means that one can consider the potential ( 137|) as a result of 
Darboux transformations applied to both 

6 2 ^ + !) 

V = x b - vx 1 + 



X 2 



and its dual (l3Uj) . In the second case one should replace the subset I C 
{1,2, ...,|M|} by its complement I* = {1, 2, ..., \M\} \ I (and thus m by 
\M\ -to). 

Proposition 5.2 The results of all possible Darboux transformations applied 
to the operator UiOi) form one of the 4 sets D±± described above. These sets 
consist o/2' A/ ' monodromy-free potentials each and are related by the dualities 
$EB\) , fijflj)- The corresponding Schrodinger equations have \M\ quasi-rational 
solutions of the form 

ipi(x) = Ri(x) exp (±x 4 /4) 

with some rational functions Ri(x). 
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Example. We saw above that the operator 

L = -D 2 + x % - 7x 2 
has two quasi-polynomial eigenfunctions 

^1,2 = {±V2x 2 + 1) exp(-x 4 /4), 

corresponding to Ai = —2\f2 and A 2 = 2\p2 respectively. 

The Darboux transformation corresponding to ip\ gives a Schrodinger 
operator with potential 

Vl = x * _ 7x 2 _ 2D 2 ln ^ = x 6 _ x 2 + «•'■' - l\/2 



(v^ + l) 2 ' 

Note that because V>i is the ground state of the operator L, the result is 
non-singular on the real line. The corresponding Schrodinger equation 

, n2 6 2 8x 2 -AV2 , A . , 

( - D+I - I+ (Vl^Tlj^ = v 

has two quasi-rational solutions: 

^ = * rl = Vl?TT exp(:c4/4) 

with A = Ai = —2\f2 and 

4^ 



2 = (£> - £>ln^i)^2 = -=V— exp(-x 4 /4) 



with A = Ai = — 2v2- The first function grows at infinity, and thus has no 
real spectral interpretation, while the second one is a genuine eigenfunction. 
Because it has a zero at x = it is the first eigenfunction above the ground 
state of the transformed potential. Note that the new ground state therefore 
is not quasi-rational. 

Similarly the Darboux transformation corresponding to ip2 gives the po- 
tential 

Applying two successive Darboux transformations (which commute) to 
both ipi and ip 2 , we come to the potential 

d 2 d 2 

Vi )2 = x 6 -7x 2 -2— ln(W(^i, V2)) = x 6 +5x 2 -2— ln(W(v / 2a; 2 +l, >/2x 2 -l)) 

16 



d 2 2 

= x 6 + 5a; 2 - 2— ln(4v / 2x) = x 6 + 5x 2 + — , 
dx z x 2 

which is of course dual to the initial potential V = x 6 — 7x 2 . 
More examples led us to the following 

Conjecture. The polynomials Pi{x) always have simple zeroes, so all 
the corresponding potentials Vj have the form {Hfy . 

We should note that there are monodromy-free potentials which are not 
of this form, as the following example shows. 
Consider a potential of the form 

V(x) = x 6 - ^x 2 + — *— + — ^— . (42) 

(x — a) 2 (x + a) 2 

The trivial monodromy conditions give two relations 

6a 5 - 2ua r = 0, 10a 3 = 0. 

2a 3 ' 8a 5 

The solution exists only when 

51 



4V15 

in which case 

8 3 
a 8 = — . 

80 
Note that the corresponding v is irrational; therefore such a potential 
can not be a result of Darboux transformations discussed above, It can not 
even have quasi-rational solutions, as will follow from the results of the next 
section. 

6 Quasi-rational solutions and Stieltjes rela- 
tions 

In this section we are going to discuss the relation between two properties: 

trivial monodromy and the existence of quasi-rational eigenf unctions. 

We start with a more general situation. Let V be a rational potential of 

the form 

N 2 

^(x) = P(s) + £-___ (43) 

i=\ ^ ' 
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^(x) = ^j -expQ(x), (44) 



where P(x) is a polynomial and Xi are distinct. Assume that the correspond- 
ing Schrodinger equation 

(-D 2 + V(x))ip = Xtp 

has a quasi-rational solution of the form 

Uf=i( x -Vj ) 

iii=l\ X ~ X i 

where Q(x) is another polynomial. The corresponding logarithmic derivative 

K N 

/ 4 rp ni . / -> rp rn . 

j=l yj i=\ l 

must satisfy the Riccati relation 

f + f 2 = V-X. 

In particular, the residues at all the poles in the left hand side must be zero. 
This gives us the following relations between poles Xi and zeroes y 3 - of any 
such eigenf unction 

K 1 - 1 

j=l Xk Vj i+k Xk Xi 

K N 

E E + Q'( yi ) = 0,l = l,...,K. (46) 

We will call these the Stieltjes relations because in the simplest case, when 
Q is quadratic and K = 0, they reduce to the classical Stieltjes relations for 
the zeroes of Hermite polynomials (see more on this in [B] ) . 

It turns out that these relations are strong enough to imply the trivial 
monodromy property for the corresponding potential V(x). 

Proposition 6.1 If the Schrodinger operator with potential ^J3fj has a quasi- 
rational eigenfunction of the form ftJ4\ ) then this operator is monodromy- free. 

The proof is elementary (cf. Lemma 2 from [B]). Let 

/ = h C + Ci(x - Xi) + c 2 {x - Xi) 2 + ... 
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be the Laurent expansion of the logarithmic derivative at a pole x = X{. The 
z-th Stieltjes relation (I4"51) means that c = 0. We have then 



f'+f 2 = ( \-c 1 (x-x i )+c 2 (x-x i ) 2 +. ..)'+( hci(x-x i )+c 2 (a;-a; i ) 2 +. 

2 2 
= t ^ - ci + (2c 2 - 2c 2 )(x - a:,) + ... = -. -^ - c t + 0(x - Xi) 2 . 

(X-Xi) 2 (X-Xi) 2 

Since /' + f 2 = V — A, we see that the trivial monodromy condition V\ — 
is satisfied. 

Corollary 6.2 The Stieltjes relations imply the locus equations for the cor- 
responding potential V(x): 

- 4 

-E fx —* S + P'{x k ) = 0,k = l t ...,N. (47) 

We will give a direct proof of this in the next section where we consider 
the time-dependent version of the Stieltjes relations. 

Note that the corresponding polynomial P(x) is defined here simply as 
the polynomial part of the rational function 



Q"(x) 



(QT+2Q'(i:— — e— ) 

\j^x- yj j^X-XiJ 



We can be more explicit in the sextic case when Q(x) = ±x 4 /4. 

Proposition 6.3 Suppose that the Schrodinger operator with potential fiJB) 
has a quasi-rational eigenfunction of the form 

^) = ^ 1 . (X " yj ' ) expeg 4 /4, (48) 

where e = ±1, then polynomial P(x) has the form 

K N 

P(x) = x 6 - e{2N - 2K - 3)x 2 + 2e(J^ y s - ^ x t )x + const (49) 

3 « 

In particular the coefficient of the quadratic term 

v = ±(2N - 2K - 3) 
is an odd integer. 
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Indeed, in that case, 



I ^_^^ rp /} / . Z—^^ up Of* . / 1 J <f rpZi ^_^^ rp rpZi 

\ iXj LJ n JU JU q I \ JU JU JU JU 



K N K N 



Note the existence of a linear term 

K N 



2e (^2vj-^2 x i) x ^ 



which does appear already in the simplest case K — N — 1, as the following 
example shows. 

Example. Consider the quasi-rational eigenf unctions of the form 



The Stieltjes relations 

1 



ip(x) — exp — x 4 /A (50) 



4 = o, 1/1 = 



xi -yi yi - xi 

imply that xf — yf, so y± = ujxi, where 

-li-v/3 

is a root of the equation u 2 + uj + 1 =0. The pole xi satisfies the equation 

1 



x 



1_ l-UJ 

and the corresponding potential is 

V(x)=x 6 -3x 2 -2(tu-l)x 1 x. (51) 

The linear term of course disappears in the symmetric case, for then 

M N 



5^ = 5Z^ = °- 
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We are going to restrict ourselves to the symmetric case from now on, allowing 
the additional ^ \ ' term at zero: 

ar 

V(x) = x 6 - vx 2 + l -^±^ + V . 2 .. . (52) 

X z ^— ' [X — Xi) 

i=l v ' 

The existence of quasi-rational solutions of the form 

tK 



,U j= i( x -yj) 

lli=l\ X ~ X i 



^(x) = x^^T 1 , r ex P ( ±a;4 / 4 ) 5 (53) 



with ji — — I or n — I + 1, leads to the following Stieltjes relations: 

K 1 ^ 1 

]T^— -]T^— ±x 3 fc + ^ = 0, fc = l,...,iV (54) 

■ =1 Xk Vj i _ Lk %k X % Xk 

K N 1 

5^_!_-^—!_± y f + ^ = I / = l,...,M. (55) 

Exactly as before, they imply the locus relations 

A 2 1(1 + 1) r 

3^i J l 

with 

v = ±(2N -2K-2/JL-3). 

Note that this is true for all values of I (and thus //) , not necessarily integers. 

The question is whether all solutions of the locus system can be found 
in such a way. In other words: is it true that any even monodromy-free 
potential of the form ( l52l) must have a quasi-rational eigenfunction? 

The answer is 'No' as the simplest example of the potential 

V{x) = x 6 - (2/ + l)x 2 + ^±^ 

x 2 

shows (see section 3 above). A more interesting example is given by 
U = x 6 - x 2 + 2/(x - a) 2 + 2/(x + a) 2 . 

The locus relation is 

12a 8 - 4a 4 - 1 = 0, 
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which gives us two possibilities: either a 4 = 1/2 or a 4 = —1/6. Since v = 
±(4 -2K-2fi- 3) = ±(1 - 2K - 2fi) = 1 and \i = or \i = 1 one can 
see that K must be either zero or 1 and we have two possible quasi-rational 

eigenf unctions: 

x — b 



^ =; — — ex P(~ x / 4 ) 

and 



x 2 — a 2 



ip= — -exp(x 4 /4). 



x 2 — a 2 



In the first case the Stieltjes relations are 



1 1 „ 11, 

a = °> r + 7T- + a 



a — b 2a —a — b 2a 

implying b = and a 4 = 1/2. In the second case we have just 

+ o 3 = 

2a 

which give the same a 4 = 1/2. Thus we see that the solution a 4 = —1/6 of the 
locus conditions is inconsistent with the Stieltjes relations, and thus the cor- 
responding monodromy-free potential has no quasi-rational eigenf unctions. 

We note that the potentials with a 4 = 1/2 are exactly those given by 
the Darboux transformations applied to V = x 6 — 7x 2 (see example in the 
section 3 above). This raised the following natural question. Suppose that 
the monodromy-free potential ( |52i) has a quasi-rational function of the form 
( |53l) . Is it true that it is a result of Darboux transformation applied to the 
operators ( |20|) ? 

We will show that this answer is negative as well. This will follow from the 
results of the next section, where we investigate the corresponding Darboux 
set in more detail. 

7 Geometry and arithmetic of the Darboux 
set 

By the Darboux set D, we will mean all possible results of Darboux trans- 
formations applied to 

V = x s _ vx 2 + l l±R. (57) 

x l 

It is the union of the sets D±± described in section 3. We will allow Iq 
(and the corresponding jjl in ( )44l) ) to be arbitrary (not necessarily integer) 
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parameters in this section. The corresponding Schrodinger operators will 
have some monodromy at x = in that case but will be monodromy-free 
elsewhere. This does not affect the relation between the locus and Stieltjes 
algebraic systems shown in the previous section. 

We know that the the operator with the potential ( [571) has a quasi-rational 
solution (1441) only if some of the numbers 



M 



±v±(2l + l) 
4 



are integers; in that case, as shown in section 3, the number of such solutions 
is equal to \M\. 

It is convenient to represent these cases on the (I, z/)-plane as the grid T 
of the lines 

v ± (21 + 1) = 4M ± 

with integer M. 




v 



M_ 



M, 



Strictly speaking, when M± = we do not have quasi-rational solutions 
but we want anyway to include the corresponding two lines (drawn dashed) 



i/±(2/ + i; 



0. 
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which we call special. The number of quasi-rational solutions is proportional 
to the distance from these lines. 

Note also that all the intersection points of the grid have I coordinates 
which are half-integers. In these cases we can have all four different types of 
quasi-rational solutions simultaneously, while at all other points in the grid 
we have only two of them. 

Now let us consider the Darboux descendants of these potentials. They 
have the form 

Vl = x e - vx 2 + ^^ - 2D 2 In P 2d (x) 
x z 

with the parameters v and I different from the original form f l57l ). 

Our first claim that these new parameters lie on the same grid T. One 
can easily check this in all the cases D±±. For example the Darboux trans- 
formations in -D+_ case, change 

v — > v — 6m, I —> I + m, 

so 

v-(2l + l) v-(2l + l) 

M = v ; -+ v ; - 2m = M - 2m, 

4 4 

where m is an integer between 1 and \M\. So we see that M remain integers, 

but with modulus \M\ decreasing, we moved on the grid closer to the special 

lines. One can check that the same change 

M -> M -2m (58) 

is true for all other types D±± as well. 

Let us look now at the degree of the corresponding polynomials P2d{x), 
which is the same as the number N of the poles in the corresponding poten- 
tial. We know that in all the cases it has the product form 

N = 2d = 2m(\M \ -m), (59) 

where 

= ±^,±(2t, + l) 



corresponds to the initial operator 

\r 6 2 , ^o + 1) / fin N 

V = x — vqx -\ . (60) 

x 1 

This means that the arithmetic of iV plays a substantial role here. In 
particular if the number of poles is 2p where p is prime then this implies 
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that corresponding m must be either 1 or p = \M \ — 1 (which is actually a 
dual case to m — 1, see section 3). So we can claim that such a potential 
must be a result of just one Darboux transformation and the corresponding 
parameters (I, v) must lie on the lines 

i/±(2/ + l) = 4(p-l). (61) 

At a generic point on these lines the initial parameters u , X are determined 

uniquely, so we have exactly p + 1 such potentials depending on the choice 

of one of the corresponding \M \ = p + 1 quasi-polynomial eigenfunctions. 

If d is not prime then the number of Darboux steps m could be any divisor 

of d = mk, where k can be assumed larger than m. In that case \Mq\ = k + m 

and we have 

{m + k)\ 

m\k\ 
potentials on the lines 

u±(2l+ 1) =4(fc-m). (62) 

Let us come back now to the question whether any monodromy-free op- 
erator with a quasi-rational eigenfunction belongs to the Darboux set or not. 
The following simple example shows that the answer in general is negative. 

Example. We consider the potentials of the form 

U = x 6 - 3x 2 + 2/(x - a) 2 + 2/(x + a) 2 
and check whether they have quasi-rational functions of the form 

x — b 

^(x) = ~o 5" ex P (~ x /4) 

x z — a z 

with a and b different from zero. The Stieltjes relations here, 
1 26 1 2a 



0, -— + „ ~ -a 3 



2b b 2 -a 2 ' 2a a 2 - b 2 

are equivalent to 

a 4 - 6 4 = 1, a 4 + 4aV + b A = 0. 

This gives 8 solutions 

4 1,2 2 

a = -, b = ea , 

1 — e z 

where e satisfies the quadratic equation e 2 + 4e + 1 = 0, giving: 

e = -2 ± V3. 

Now we claim that the corresponding potentials cannot be the result of 
Darboux transformations. Indeed in that case d — 1, som must be 1 and 
\Mq\ = 2. On the other hand, since I = 0, v = 3 we have M = ^p = 1, 
which must be equal to Mq — 2m = Mq — 2, which is clearly impossible. 
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8 Time-dependent Stieltjes relations and spe- 
cial solutions of the Calogero-Moser prob- 
lem 

We now consider the time- dependent Schrodinger equation 

If V(x, t) = x 6 — vx 2 + - 2 , then this has a linear space of dimension M of 
quasi-rational solutions, with exponential time-dependence: if ip l satisfies 

V{x)W = \^\ 



dx 2 



then ip c = Ylj=i c j^ exp(— zAjt/2) evidently solves fl63l) . This solution, which 
has simple but non-trivial time-dependence, may be used, as above, to con- 
struct a Darboux transformation 

^ 1 ,...,c M =^-2^1n(W^ 1 (<)). 

The moving poles of the new potential must thus satisfy the monodromy 
conditions; in particular they must solve the corresponding Calogero-Moser 
system. 

As an example let us consider the N = 2 case. Putting exponential 
time-dependence into the eigenfunctions discussed above, we construct a 
2-dimensional space of quasi-polynomial solutions of the time-dependent 
Schrodinger equation: 

namely 

^ c = ( Cl (V2x 2 + l)e iV2t + c 2 (V2x 2 - l)e~ iV2t ) exp(-a; 4 /4). 
The Darboux transformation associated with this leads to a new potential 

VW=x 6 -x 2 -2^ln^ c ) 

6 2 , , X 2 +X 2 (t) 

= x° - x 2 + 4- 



{x 2 -x 2 (t)y 
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where 



isftt _ -iV2t 

X\t) = 



—l- 



^/2( Cie iV2t + C2 e-^t) 
tan( v / 2(t-t )) 

We may easily check in this simple case that the points 

Xl (t)=X(t), 

x 2 (t) = -X(t) 
indeed satisfy the Calogero-Moser system 

E Q 5 

[Xi-Xjf 

From the locus condition, such a result in fact holds for all potentials of this 
kind, as will be shown more directly below, where for simplicity we assume 
the parameter / = 0. 

As in the time-independent case, Darboux transformations take a quasi- 
polynomial solution of the time-dependent Schrodinger equation with an even 
sextic potential 

V^\x) = x 6 - (2ra + 3)x 2 

into a quasi-rational solution 

wm+l/ /(OK 
dm) = W j=l ^ C J j 

rci,...,c m ;c m+ i //(OK' 

which satisfies the Schrodinger equation with the rational potential 

K£lc m =^~ (2n + 3):r 2 - 2^ ln(Wf =1 (^f) 

<9 2 
= a; 6 - (2ra + 3 - 6m)x 2 - 2— ln(W™ 1 (0| ; o) ). 

Here, as before, the 0c- are the rational factors of the corresponding quasi- 
rational solutions ^ Cj . 

As before, this solution may be described in terms of its poles and zeroes 

as 

^ m) = gj=lfe - *j)°' exp( _ x 4 /4) exp(i/(t))> 
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where the multiplicities aj and f3k are all positive integers. Here the points 
Xk, Xk are functions of time, while the time-dependent phase /(£) is not yet 
determined. We identify the numerator and denominator of the rational part 
with the polynomial factors of the two Wronskians. 

z 
W m+1 = l[(x - Xj p = WfftfW), 

i=i 

p 
W m = l[(x - xtf* = Wf =1 (0W). 
fe=i 

We note that the numerator and denominator may have common factors. 
The Schrodinger equation is now 

. d W m+1 - W m+1 ^ 
% dt W m ; W m 

.i ( | _ ,3 f + i x « _ {N _ 3m + mx > _ JL ln(Hy ) i^±i, 

and the Hirota form is hence: 

.dW m+1 dW m ■ 

i gt W m - iW m+1 — fW m+1 W m 

l ( d 2 W m+Ujr n dW m+1 dW m , TT7 a 2 W, 



+-( ' " W m - 2 '"^ ' '" + W, 



m n 



m+1" 



2 ax ax ax oar 

3 9WWi aW m , 

-x ( — r W m - W m+1 - 



OX OX 

= x 2 {N-3m)W m+1 W m . 
Now, on dividing the whole equation by W m W m+ i, this becomes: 

z P ' 

ECXjXj .y-^ (3 k X k x 

— L ^ + i }_^ —~ f 
ry ry . ^^^^4 ry ry 

j=l J fc=l 



2 I -^ 2-^ (x — Xh)(x - Xj,) 2^2^ ( x - xAtx - Xk) ^ ^ 



1 fc= l )fc? 4j ( X X i)( X X fe) i=1 fc=1 ( X X 3>\ X X k) j=lk=hk ^j( X X J)( X X k) 

+x 3 ( y- _^_ _ y _A_) = x 2 (Ar _ 3m) 



j=l J fe=l 
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Matching coefficients of x 2 , we again get ^2 j=1 olj — Sj=i Pj = N — 3m. The 
coefficients of x must all vanish, by parity, and the coefficient of x° can now 
be balanced by requiring 

z p 

-/ +x;«^ -£&**= °- 

i=i i=i 

It remains to check the residues at the zeroes Xj and poles 5;^ which give, 
respectively: 

Z P a 

- lajXj + 2_^ 7 — ^ ~~ /_, 7 — - \ + x j a j = °> (64) 

and 

i/te- £ 7^^ + E 7^T + ^ = 0. (65) 

• 1 _L1 \%k %j) . \Xk %j) 

These are the dynamical Stieltjes relations we seek. 

To simplify the notation, if we look at the union of the sets of poles and 
zeroes, points Zi, each with a positive or negative exponent 7$, 

Zi Xj, 7« ®-ii ^ ij ■ ■ ■ j "j 

^ = 5j_ z , ji = -pi_ z , i = Z + 1,...,Z + P, 
then this set of relations simplifies: 

z+p 



k=l,kft ^ k) 



It is instructive to differentiate these relations with respect to time. We find: 

z+p 



Zj =1 



k=i,k?j { J k> 
On substituting (1661) into this, we obtain 

z+p 

Zi = i 

k=l,k=£j 



1 " ^A*j-*h) 



Z+P Z+P 

Ell . 3 . • V^ ll . • 3\ 
--iz.+i y - + iz k ) 
(zj -zi) 3 ^ (zk - zi) 

Z+P 

-z%z){-% J2 {7 lk _ 7 ^ -*4)- 

l=l,l& Kh l) 
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This may be rearranged as 

z+p 



y —> (Zj - z k f 



k=l,k^j v ° 

Z+P Z+P 



z k> 



(Zj — Zh) 2 [z,\ — Zi) 3 (Zh — Zi) 

Z+P 



-M( E 



--IZj). 



(Zj — Zi) 

This further simplifies, after some elementary manipulation, to 

z+p , s z+p z+p z+p 

z 3~ 1^ ( z .- z )3 — 6z i~ 1^ 1^ ( z .- z )(z--z)(z -z)~ 2-" 7k{lZj+Zk) 

k=l,kft ^ j k) k=l,k^j l=l,l^j,k ^ j fe ^ j l '\ k l > k=l,k^j 

Now in this long expression the double sum 

z+p z+p 

Ikli 



E E 



is identically zero. The final term 

z+p 
^2 lk(2zj + z k ) 

k=l,k^j 

is equal to 



z jK 



Z+P Z+P 

,(2 ^2 Ik- 37j) + Y^ 7fc*fc)- 
fc=i fc=i 

In the special case where all zeroes and poles are simple, 7^ = ±1, and 
X)fc=i 7fc = ^ — -P) so this term becomes Zj(2(Z — P) — 3"fj). Further, from 
the quadratic term above we saw Z — P = N — 3m. The term Ylk=i IkZk is 
supposed to be identically zero by parity. 

Thus, remarkably, the system decouples, giving: 

Z 2 

x s = £ — ——-:**- Xj (2(Z-P)- 3). (67) 

k=\,k+3 { j k) 

for the zeroes, and 

z+p 

k=l+Z,k^j ^ j k > 

for the poles. 
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Proposition 8.1 The Stieltjes relations (64\ [6d*\) withai = /3j = 1 imply that 
the zeroes and poles of solutions of the non- stationary Schrodinger equation 
satisfy the uncoupled Calogero-Moser systems in different sextic potentials 
with the Hamiltonians 

i=l i=l i=l i=l j=l,j=£i v ]J 

for the zeroes, and 

i=l i=l i=l i=l j=l,j& v ■ J/ 

/or i/ie poles. 

Thus the first order equations ( 1661) again imply that the simple poles and 
zeroes satisfy separate Calogero-Moser systems, in even sextic potentials, but 
with different 0(x 2 ) terms. We may thus regard (1661) as being a Backlund 
transformation for Calogero-Moser systems of this type. This does not imply 
the system is integrable however, for the Backlund transformation does not 
depend on an arbitrary constant. Further, it can only be iterated finitely 
many times. 

Alternatively, the system of dynamical Stieltjes relations may be under- 
stood as a canonical transformation. Taking all the poles and zeroes to be 
simple, if the momentum variable conjugate to the zero Xj is pj, and that 
conjugate to Xj is pj, the relations read: 

Z 1 P 1 

Pi = - 1 E 7^T^ +l E^T^- l 4 (69) 

k=i,k?j { 3 k) fe=i l 3 k) 

and 

Z 1 P 1 

Pk = i'^2j^ r~ i ^2 ~F ^ ~~ ix3 k = °- ( 70 ) 

j=i \ Xk — x i> j=\.j+k ^ Xk ~ X °> 



If we define 

lj=l \.lk=l,k>j\ jj 3 — -^k) ll m =i lln=l,n>m 



. , , 1 lj=l 1 lfc=l,jfc>j \ X j X k) H m =l [[ n =l,n>m,\ X m X n 

S = —ilog 



U j= i U k =i,k^j( x i - x k) 



Z 4 P ~4 

i=i fc=i 
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then the Stieltjes relations (I6TH 1701) become: 

»• = If- < 72 > 

* ^ "f • (73) 

They are thus seen to give a canonical transformation with generating func- 
tion S. Note that dropping i in these formulae corresponds to a change of 
sign in the potential and thus to an attractive, rather than repulsive, pairwise 
potential. 

This transformation is straightforward when the numbers Z and P of 
zeroes and poles are equal, but otherwise we may still understand it as gen- 
erating a mapping from simple solutions of a Calogero-Moser system, with 
few degrees of freedom, to more complicated solutions, with more degrees 
of freedom - of course such solutions will be special, with some algebraic re- 
lations between the new coordinates and momenta. The rational potentials 
constructed above correspond to the most special ones of all these, which are 
connected by a finite chain of such canonical transformations to the 'empty' 
Calogero-Moser system with no particles at all. Since it is elementary to show 
that this transformation preserves the Hamiltonian, H = H, all the solutions 
of Calogero-Moser systems constructed here, related by finitely many such 
transformations to one with no degrees of freedom, must have zero energy. 
We believe that the same should be true for other Inozemtsev integrals as 
well. 

One immediate consequence of the fact that the dynamical Stieltjes re- 
lations imply that the simple poles and zeroes separately satisfy Calogero- 
Moser equations is that the analogous stationary Stieltjes relations must im- 
ply that the poles and zeroes satisfy stationary Calogero-Moser equations in 
the appropriate potentials; these are precisely the locus conditions for these 
potentials. 

9 Concluding remarks 

We have seen that the monodromy-free potentials with sextic growth are 
described by the algebraic locus equations and are related to equilibria of 
the corresponding version of the Calogero-Moser problem, which is known to 
be integrable [S]. So, is the locus system (157)1) 'integrable' ? In which sense ? 
We have shown that Darboux transformations are not enough to describe 
the whole locus, even if we assume the existence of quasi-rational solutions. 
The geometry and arithmetic of the corresponding potentials seems to be 
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quite interesting as well as the algebraic geometry of the corresponding (zero- 
dimensional) affine varieties. In particular, the Stieltjes relations can be 
considered as correspondences (in the sense of algebraic geometry) between 
two different locus varieties, which is a partially defined multivalued algebraic 
map. Although they can be considered as natural analogues of the Backlund 
transformations their existence is not enough for 'integrability' and holds for 
any growth at infinity (see section 6 above). 

But the general structure of the sextic locus still remain largely unclear. 
In relation with Inozemtsev's results [8] we note that the general solution of 
the Calogero-Moser system in a sextic potential is given in terms of some 
matrix analogues of the elliptic functions and are expressed in terms of the 
high genus Riemann ^-functions. This means that an elementary description 
of all corresponding equilibria seems to be unlikely in agreement with our 
result. 

In that sense the situation looks similar to the case of the elliptic poten- 
tials studied by Gesztesy and Weikard [21], who showed that all the corre- 
sponding potentials must be finite-gap in the sense of [22]- The problem of 
finding an effective description of the elliptic finite-gap potentials was first 
raised by S.P. Novikov in the 70s, who was inspired by the first examples 
found in his work with Dubrovin [23]. This area was revitalised after the 
appearance of the famous Trebich-Verdier preprint [2lj , which showed that 
there are examples, which can not be found using the KdV flows from the 
classical Lame potentials 

V(x) = m{m + l)p(x), 

where p(x) is the Weierstrass elliptic function. In spite of the numerous ef- 
forts in this direction (see e.g. Enolskii-Eilbeck [25], Taimanov [26], Smirnov 
|27j ) the problem seems to be far from a satisfactory solution. 

Our sextic problem can be considered as a simple model for this important 
problem and thus deserves further investigation. It could also very interest- 
ing from the analytic-differential point of view, see the recent paper [2H] by 
Eremenko, Gabrielov and Shapiro. 
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